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Abstract 

In cosmology an important role is played by homogeneous and isotropic 
solutions of the Einstein-Euler equations and linearized perturbations of 
these. This paper proves results on the asymptotic behaviour of scalar 
perturbations both in the approach to the initial singularity of the back- 
ground model and at late times. The main equation of interest is a linear 
hyperbolic equation whose coefficients depend only on time. Expansions 
for the solutions are obtained in both asymptotic regimes. In both cases 
it is shown how general solutions with a linear equation of state can be 
parametrized by certain functions which are coefficients in the asymptotic 
expansion. For some nonlinear equations of state it is found that the 
late-time asymptotic behaviour is qualitatively different from that in the 
linear case. 



1 Introduction 

Astronomical observations allow information to be collected about the distribu- 
tion of matter in the universe. This distribution contains structures on many 
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different scales. Astrophysicists would like to provide a theoretical account of 
how these structures formed. In particular, cosmologists would like to do this for 
structures on the largest scales which can be observed. This means for instance 
giving an explanation of the way in which galaxies cluster. The most power- 
ful influence on the dynamics of the matter distribution on very large scales is 
gravity. The most appropriate description of the gravitational field in this con- 
text is given by the Einstein equations of general relativity. It is also necessary 
to choose a model of the matter which generates the gravitational field. A fre- 
quent choice for this is a perfect fluid satisfying the Euler equations. Thus, from 
a mathematical point of view, the basic object of study is the Einstein-Euler 
system describing the evolution of a self-gravitating fluid. This is a system of 
quasilinear hyperbolic equations. 

The standard cosmological models are the Friedmann-Lemaitre-Robertson- 
Walker (FLRW) models which are homogeneous and isotropic. This means in 
particular that the unknowns in the Einstein-Euler system depend only on time 
and the partial differential equations reduce to ordinary differential equations. 
With appropriate assumptions on the fluid these ODE's can be solved explicitly 
or, at least, the qualitative behaviour of their solutions can be determined in 
great detail. When it comes to the study of inhomogeneous structures, however, 
the FLRW models are by definition not sufficient. Since fully inhomogeneous so- 
lutions of the Einstein-Euler system are difficult to understand a typical strategy 
is to linearize the system about a background FLRW model. Under favourable 
conditions the linearized perturbations could give information about the evolu- 
tion under the Einstein-Euler system of initial data which are small but finite 
perturbations of those for the FLRW background. 

Linearization about a highly symmetric solution is a classical practise in 
applied mathematics. For some examples see [H], |5] and [5]. It should be 
noted, however, that there is an unusual feature in the case of the Einstein-Euler 
system which has to do with the fact that these equations are invariant under 
diffeomorphisms. This is related to the fact that the only thing that is of physical 
significance are equivalence classes of solutions under diffeomorphisms. Since 
it is not known how to develop PDE theory in a manifestly diffeomorphism- 
invariant way this leads to difficulties. There is a corresponding equivalence 
relation on linearized solutions. Different linearized solutions are related by the 
linearizations of one-parameter families of diffeomorphisms, which are known 
in the literature on cosmology as gauge transformations. In the end what is 
interesting is not the vector space of solutions of the linearized equations but 
its quotient by gauge transformations. It is useful to represent this quotient 
space by a subspace. This is what is known in the literature on cosmology as 
gauge-invariant perturbation theory. This subject would no doubt benefit from 
closer mathematical scrutiny but that task will not be attempted in the present 
paper. 

Instead the following pragmatic approach will be adopted: take an equation 
from the astrophysical literature on cosmological perturbation theory and anal- 
yse the properties of its solutions. As a basic source the book of Mukhanov [8] 
will be used. The notation in the following will generally agree with that of 
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[S]. It is standard to classify cosmological perturbations into scalar, vector and 
tensor perturbations. These terms will not be defined here. It should be noted 
that scalar perturbations play a central role in the analysis of structure forma- 
tion. This motivates the fact that the results of this paper are concerned with 
that case. After a suitable gauge choice scalar perturbations are described by 
solutions of a scalar wave equation for a function <I> which corresponds, roughly 
speaking, to the Newtonian gravitational potential. In order to get definite ex- 
pressions for the Einstein-Euler system and its linearization about an FLRW 
model it is necessary to choose an equation of state p = /(e) for the fluid. Here 
e is the energy density of the fluid and p its pressure. A case which is partic- 
ularly simple analytically is that of a linear equation of state p = we where w 
is a constant. For physical reasons w is chosen to belong to the interval [0, 1]. 
In fact the condition w > is necessary in order to make the Euler equations 
hyperbolic. The case w — 0, known as dust, is somewhat exceptional and does 
not always fit well with the general arguments in the sequel. Since, however, 
dust frequently comes up in the literature on cosmology it is important to in- 
clude it. In those cases where the general argument fails for dust this will be 
pointed out. 

For a linear equation of state as just described the equation for <& is 

^ 6(1 + ^^1^,^ (1) 

Here a prime stands for The time coordinate 77 belongs to the interval (0, 00). 
The spatial variables, which will be denoted collectively by x, are supposed to 
belong to the torus T^. Thus periodic boundary conditions are imposed. The 
Laplacian is that of a fixed flat metric on the torus. Its expression in adapted 
coordinates agrees with that for the usual Laplacian on R^. As a consequence 
of standard theory for linear hyperbolic equations this equation has a unique 
solution on the whole time interval (0, cxd) for appropriate initial data given at 
a fixed time t] = rjo > 0. These are the restrictions of $ and to r] = rjQ. 

In the following, after some background and notation has been collected in 
Sect. [21 the asymptotics of solutions of equation ^ is studied in the regimes 
ri ^ and rj ^ 00. Theorems and proofs for the first of these cases are given 
in Sect. [3] (Theorems [T] and [2]) and for the second in Sect. H] (Theorem S]) . It 
is shown how all solutions can be parametrized by asymptotic data in either of 
these regimes. These are alternatives to the usual parametrization of solutions 
by Cauchy data. An interesting feature of the expanding direction 77 — > cx) is 
that the main part of the asymptotic data is a solution of the fiat space wave 
equation W" = wAW. Many of these results can be extended to more general 
equations of state. This is the subject of Theorem [3] of Sect. [3] (limit 77 ^ 0) 
and Sect. 5. It is found that for equations of state with power law behaviour 
p ~ e^^*^ at low density there is a bifurcation with a fundamental change in the 
asymptotic behaviour a.t a — ^. 
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2 Preliminaries 



As outlined above, we study perturbations of FLRW cosmological models which 
are spatially flat and have spatial topology. The spacetime being perturbed, 
which we refer to as the background, is described by a metric of the form 



a' 



{-djf + da;2) (2) 



on (0,00) X . Here dx^ indicates the flat metric on and the scale factor 
a = a{ri) is a non-decreasing function of the conformal time 77. We use x to 
indicate points on T^. The signature used here is the opposite of that used by 
Mukhanov [S] but all the equations required in the following are unaffected by 
this change. 

We make use of the perfect fluid matter model, described by the pressure 
p and energy density e of the fluid. In order to specify the matter model com- 
pletely, one must provide an equation of state p = f{e). Under this ansatz, the 
Einstein-Euler equations reduce to a coupled system of ODEs for a and e: 

a" = ^(e-3/(e))a3 (3) 
e' = -m{e + f{e)). (4) 

As mentioned in the introduction, ( )' indicates a derivative with respect to 
77. Here G is Newton's gravitational constant and Ti is the conformal Hubble 
parameter, given by 7i = a~^a'. We note the following useful relation (known 
as the Hamiltonian constraint) 

„ ,9 8itG 9 

= —ah. (5) 

For a linear equation of state /(e) = we, solutions a{ri) of ^ are explicitly given 

by 



, N /X 2/(l-|-3t«) 

a[r]j I V 



(6) 



for some arbitrarily fixed 770 G (0,oo). As the scale factor a vanishes as 77 ^ 0, 
the spacetime develops a curvature singularity in that limit, which is known as 
a "big-bang" type singularity and is viewed as being in the past of 7^9 . Likewise 
the limit as ry ^ cx) is referred to as "late times" as it corresponds to the distant 
future of 770. Note that spacetimes described by these models are expanding, in 
the sense that the scale factor is an increasing function of 77. Note also that, since 
e' is negative, large values of ry correspond to small values of e and vice-versa. 

We study behavior near the singularity and at late times for those pertur- 
bations to the metric ([2]) which are of the type usually referred to as scalar 
perturbations. They satisfy evolution equations obtained by linearizing the 
Einstein equations about the FLRW background. For the perfect fluid matter 
model all such perturbations can be described, up to gauge freedom, by a single 
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function ^{t],x). Using a certain gauge, the conformal-Newtonian gauge, the 
metric takes the form 



a 



(7) 



up to an error which is quadratic in the expansion parameter A. The first order 
perturbation satisfies the linearized Einstein-Euler equations provided 



$" + 3 (f + /'(e)) + 3 ( /'(e) - ^ ) - /'(e) A$ = 0, (8) 



where A is the Laplacian for the flat metric on T^. For a derivation of this 
equation we refer the reader to §7.2 of [8]. The corresponding perturbations to 
the energy density, denoted by Se, are determined by 



and thus can be computed once ^ is understood. 

The quantity /'(e) represents the square of the speed of sound for the fluid. 
For physical reasons we require that /' always take values in the interval [0, 1] 
i.e., that the speed of sound be real and not exceed the speed of Hght. A special 
case of particular interest is that of a linear equation of state p = we. In this 
situation the speed of sound is constant and equation dH) reduces to (H]). Before 
the asymptotics of solutions of ([8]) can reasonably be studied a prerequisite is a 
theorem which guarantees global existence of solutions on the interval (0, oo). In 
order to get this from the standard theory of hyperbolic equations it is necessary 
to assume that /' never vanishes. In the following it is always assumed that 
this holds except in the special case of dust which is discussed separately. 

Our analysis below relies on establishing a number of energy-type estimates 
for solutions to ([8]). As the coefficients of this linear equation depend only on rj, 
any spatial derivative of $ satisfies the same equation. Thus any estimate we 
obtain for or V$ (the gradient of $ with respect to the fiat metric on 

T^) holds also for all spatial derivatives of those quantities. One may then make 
use of the Sobolev embedding theorem in order to establish pointwise estimates. 
We also make use of the Poincare estimate which implies that quantities having 
zero (spatial) mean value are controlled in by the norm of their (spatial) 
gradient. 

Each of these norms is defined on the ?7-constant "spatial" slices of (0, oo)xT^ 
with respect to the fiat (?7-independent) metric induced on by viewing as 
a quotient of Euclidean space. All integration on is done with respect to the 
corresponding volume element which we suppress in our notation. We generally 
suppress dependence of functions on the spatial variable x, except in situations 
where the inclusion of such dependence provides additional clarity. When nec- 
essary, we denote Cartesian coordinates on hy x — (x'); the corresponding 
derivatives are denoted di. 
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3 Asymptotics in the approach to the singular- 
ity 

The purpose of this section is to analyse the asymptotics of solutions of ([1]) in 
the limit — > and to give some extensions of these results to more general 
equations of state which need not be linear. Define = ^ (f+fuj) • ^o^e that v 
belongs to the interval [1,5/2]. 

Theorem 1. Let ^{rj) be a smooth solution of ^ on (0,oo) x T^. Then there 
are coefficients ^k,i with k > ~1v belonging to an increasing sequence of real 
numbers tending to infinity and I G {0, 1}, smooth functions on T^, such that the 
formal series X]fc(*^fc,o + ^k,i log?7)'7'^ is asymptotic to ^{r]) in the limit rj —> 
in the sense of uniform convergence of the function and its spatial derivatives 
of all orders. All coefficients can be expressed as linear combinations of <^-2v,0! 
4>o,o o,nd their spatial derivatives. If v is not an integer then all coefficients with 
I = 1 vanish. For any value of w the coefficients ^k,i with I — 1 and k < 
vanish. 

In more detail, ^k,o rnay only be non-zero when k is of the form —2v + 2i or 
2i for a non-negative integer i while ^k,i rnay only be non-zero for k of the form 
2i with i a non-negative integer. These coefficients are related by the following 
equations: 

k{k + 2v)^kfl = wA$fc_2,o - (2fc + 2v)^k,i (10) 

and 

fc(fc + 2i/)$fc,i = wA$fe_2,i- (11) 

Proof. The basic tool which allows the solutions to be controlled is provided by 
energy estimates. Let 

Ei{ri)^ \ I Wm^+w\VH^)\\ (12) 



It satisfies the identity 
d 

drj 



Tf{^^+^)E^(T^)\={2y + l)7^^^+^[ u;|V$(?7)|2. (13) 



Since the right hand side is manifestly non-negative it can be concluded that 
if an initial time r/o is given then jf''^'^^^^^ Eiijf) is bounded for 77 < 779. Any 
spatial derivative of $ satisfies the same equation as <i>. Thus corresponding 
bounds can be obtained for the norms of all spatial derivatives. Applying 
the Sobolev embedding theorem then provides pointwise bounds for <i> and its 
spatial derivatives of all orders in the past of a fixed Cauchy surface. These 
estimates can now be put back into the equation to obtain further information 
about the asymptotics. To do this it is convenient to write ([Ij in the form 

4 h'"+'$'(77)l = 772'^+iu;A$(77). (14) 
drj 
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It can be deduced that 



rno i-n 



rno rv 1 

, / ^2^+iA$(C)dC + w' / C^''+^A$(C)dC . (15) 
Jo Jo 



The bounds already obtained guarantee the convergence of the integrals. This 
formula allows the asymptotic expansions to be derived inductively. Using the 
fact that the second integral is 0{r]^) already gives a one-term expansion for 
and this can be integrated to give a one-term expansion for $. Analogous 
expansions can be obtained for all spatial derivatives of $ in the same way using 
the corresponding spatial derivatives of (fT5|) . When an asymptotic expansion 
with a finite number of explicit terms is substituted into the right hand side of 
an expansion for <&' (and thus by integration for <I>) with additional explicit 
terms is obtained. If the last explicit term in the input is a multiple of rjP with 
p < —2 then there is one new term in the output and it is a multiple of 
If the last explicit term is a multiple of 77^^ there is one new term and it is 
a multiple of log 77. If the last explicit term is a multiple of log 77 then there 
are two new terms, one a multiple of 77^ log 77 and one a constant. If the last 
explicit term is 77P or 77P log 77 with p > —2 then there is one new term and it is 
a multiple of 77^+^ or 77^*+^ log 77 respectively. These statements rely on the fact 
that when any of the terms in the asymptotic expansion is substituted into the 
last integral in psp the power —1 never arises. These remarks suffice to prove 
the first part of the theorem. The resulting series is by construction a formal 
series solution of the original equation. Comparing coefficients gives the rest of 
the theorem. □ 

Note that the only two values of w in the range of interest where logarithmic 
terms occur in the expansions of the theorem are w = ^ and w = 1 corresponding 
to 1/ = 2 and i' ~ 1 respectively. The two cases of most physical interest, w — 
(dust) and w — ^ (radiation), are free of logarithms. In the case w = most 
of the expansion coefficients vanish and the two non-vanishing terms define an 
explicit solution which is a linear combination of two powers of 77. 

The relative density perturbation is given by 

^ = -2$ - 2^-1$' + -7i:-2a$. (16) 
6 3 

Now Ti. = ^j^^g^^^ . Substituting this relation and the asymptotic expansion for 
$ into the expression for the density perturbation gives: 



6e X ^ 



-(fc(l + 3w) + 2)$fc,o - (1 + SH^fc.i + + 3u;)2A$fe_2,o 

6 



+ {~{k{l + 3w) + 2)*fe.i + i(l + 3w)^A<i>k-2,i logr/ 

o 



(17) 
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The relations in Theorem [T] place no restrictions on the coefficients ^-2v,o 
and $0,0 and so it is natural to ask if these can be prescribed freely. In other 
words, if two smooth functions on T'^ are given, is there a smooth solution of 
the equations in whose asymptotic expansion for r/ — s- precisely these func- 
tions occur as the coefficients $-2i/,o and $0.0? The next theorem answers this 
question in the affirmative. Since the proof is closely analogous to arguments 
which are already in the literature it will only be sketched. 

Theorem 2. Let ^1 and ^2 be smooth functions on . Then there exists a 
unique solution of ([T]) of the type considered in Theorem]^ with $_2i/,o — ^1 



and $, 



0,0 



Proof (sketch). The proof of this theorem uses Fuchsian techniques. It imple- 
ments the strategy applied in [S] to prove theorems on the existence of solutions 
of the vacuum Einstein equations belonging to the Gowdy class with prescribed 
singularity structure. In the present situation some simplifications arise in com- 
parison to the argument for Gowdy due to the fact that the equation being 
considered is linear. The procedure is to first treat the case of analytic data and 
then use the resulting analytic solutions to handle the smooth case. To reduce 
the equation to Fuchsian form the following new variables are introduced. First 
define a function v{ri,x) by the relation 



$(r/) = '^iTj 



-2v 



E 

-2iy<fe<0 



^kfiif + $0,1 logr/ + *2 + 



(18) 



Here it is assumed that the consistency relations (|10l) hold for —2v <k <Q. As 
a consequence of these relations and the original equation, v satisfies 



2v 



V 



—v' - wAv = wA^o^ilogri + wA'i^2 + w ^ A^kfiv'' (19) 



0<fe<2 



Note that the last sum will contain one non- vanishing term for h' not an integer 
and none for i' an integer. Denote the right hand side of ^9]) by Q. This 
equation can be reduced to a first order system by introducing new variables 
= rjv' and z)' = rjdiV. Let V be the vector-valued unknown with components 
{v,v^,v'^). Then the first order system is 

r]d^V + NV = ri'^firi, V, DV) 



where 



N 



-1 
2v 




and / 







ri^-<d,{v° + v) 



(20) 



(21) 



Here C is any positive real number less than one and DV denotes the col- 
lection of spatial derivatives of V. It will be shown that this equation has a 
unique solution v which converges to zero as 77 — > 0. Initially we assume that 



8 



the functions and ^2 are analytic. Then results proved in [7] can be applied. 
See also section 4 of [1] for some further information on these ideas. One of 
the hypotheses required is that / is regular in the analytic sense defined there. 
What this means is that / and all its derivatives with respect to any argument 
other than rj are real analytic for rj > and extend continuously to 77 = 0. 
The other hypothesis is that the matrix exponential cr^ should be uniformly 
bounded for all < cr < 1. This follows from the fact that is diagonalizable 
with non-negative eigenvalues. 

To extend this result to the smooth case more work is necessary. The basic 
idea is to approximate the smooth functions ^'i and ^'2 by sequences of analytic 
functions {'^i)n and (^2)n, apply the analytic existence theorem just discussed 
to get a sequence of solutions Vn of the Fuchsian system and then show that 
Vn tends to a limit 1^ as n 00. The function V is then the solution of the 
problem with smooth data. To show the convergence of Vn suitable estimates 
are required and in order to obtain these the Fuchsian equation is written in an 
alternative form which is symmetric hyperbolic. This rewriting is only possible 
for w but for w — the system, being an ODE, is already symmetric 
hyperbolic and so the extra step is not required. In general a simplification 
of the system is achieved by introducing a new time variable by i = 77^ and 
rescaling / by a factor Then the system can be written as 



tA°dtV + tA^djV 



MV = tg{t, V, DV) 



(22) 



where 



M = 



-1 
^ 











" " 




2-C 


, 9 = 


t—Q 








(23) 



and the other coefficient matrices are given by A" = diag(l, — ,/) and 



A^ = 








-r 



H—e, 




(24) 



with Cj the jth standard basis vector in . This is a symmetric hyperbolic 
system. A disadvantage is that in passing from N to M positivity is lost. 

The fact that M has a negative eigenvalue can be overcome by subtracting 
an approximate solution from v to obtain a new unknown. Expressing the 
equation in terms of the new unknown leads to a system which is similar to that 
for V but with M replaced hy M + nl for an integer n. For n sufficiently large 
this means that the replacement for M is positive definite. With this choice 
the system is both in Fuchsian form and symmetric hyperbolic. The necessary 
approximate solution can be taken to be a formal solution of sufficently high 
order as introduced in section 2 of [S]. The fact that the system is symmetric 
hyperbolic leads to energy estimates which can be used to prove the convergence 
of the sequence of analytic solutions to a solution corresponding to the smooth 
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initial data, thus completing the proof of the existence part of the theorem. 
Uniqueness can be proved using an energy estimate as has been worked out in 

M- □ 



It would presumably be possible to extend the above results to the case 
that the data are only assumed to belong to a suitable Sobolev space. An 
alternative approach to doing so would be try to apply ideas in the paper [6] of 
Kichenassamy. 

The proofs just presented have been strongly influenced by work on Gowdy 
spacetimes. For a special class of these, the polarized Gowdy spacetimes, the 
basic field equation is Pu + t~^Pt = Pxx- Evidently this is closely related to H]) 
although they are not identical for any choice of w, even if attention is restricted 
to solutions of ([T]) depending on only one space variable. The energy arguments 
above were inspired by those applied to the polarized Gowdy equation in [3]. 
The following analogue of Theorem 2 is a special case of a result in [9j . If smooth 
periodic functions k{x) and lo{x) are given with k everywhere positive there is 
a smooth solution of the polarized Gowdy equations which satisfies 



as t ^ 0. It is plausible that the positivity restriction on k, while very impor- 
tant for general (non-polarized) Gowdy spacetimes, should be irrelevant in the 
polarized case. It turns out that following the arguments used above to analyse 
(fT|) allows this intuition to be proved correct. 

One way of attempting to reduce the polarized Gowdy equation to Fuchsian 
form is to mimic and write P — k\ogt + uj + v. This fails because the 
analogue of the matrix N has v replaced by zero. Thus the matrix has all 
eigenvalues zero and includes a non-trivial Jordan block. To access the Fuchsian 
theory in the analytic case the expansion for P may be replaced by 



for a small positive 5. With this modification the reduction procedure applied to 
([D) gives a Fuchsian system. It can be concluded that k and uj can be prescribed 
in the case that they are analytic. Once this has been achieved the smooth case 
can be handled just as in the proof of Theorem 2. 

It will now be shown that some of the results which have been proved for a 
linear equation of state can be extended to more general equations of state. In 
the discussion which follows it will be convenient to exclude the case of a linear 
equation of state which has been treated already. This in particular excludes 
dust so that by our general assumptions /' never vanishes. In this case we 
consider solutions to ([5]) rather than ([T]). Choose an initial time 770 and for a 
given background solution let £(770) — Co- From the condition that /'(e) < 1 it 
follows that 



P{t,x) = k{x)\ogt + Lo{x) + o{l) 



(25) 



P = k\ogt + LO + t^v 



(26) 



1/2 



A:= sup /'(e) 



(27) 



(eo.oo) 



e 
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is strictly positive and finite. It will be assumed in addition that the equation 
of state satisfies the condition 



sup 

(co.oo) 



< oo. 



(28) 



Using the fact that A > it follows that there exists a positive number A 
satisfying the following three inequalities: 



2a|- >3(e 
ae 



/(e)) 



de2 ' 



(29) 



4A^ 



6(l + /'(e))+ 1 + 



3/(e] 



+ 6(1 + /'(£)) 1 + 



3/(e^ 



- A- 



A2 - 3 /'(e) 



fie) 



> (30) 



and 



A > 3(1 + /'(£)). (31) 

That (QSl) can be satisfied follows from The fact that /'(e) and /(e)/e are 
bounded means that the first term in the expression on the left hand side of (|30p 
dominates the other terms for A sufficiently large and so the second condition 
on A can also be satisfied. The constant A can be chosen to satisfy ([31]) since 
the right hand side of that inequality is bounded. Note for comparison that for 
a linear equation of state A = 0. In that case A can be taken to be the larger 
root of the expression obtained from the left hand side of (|30p by omitting the 
term containing A. This root is 3(1 + w). Define the following generalization of 
the energy functional (fT^ : 

E2{l) = \l |$'(,7)P + /'(e)|V<f(r])P+A7^2|cj>(^)|2. (32) 

(Note that we suppress the dependence of e and H on rj.) A computation shows 
that if a denotes the scale factor, then due to the inequalities (|29l) - ((3T|) 

^ [a^'E^iv)] > 0. (33) 

In more detail, computing the time derivative of a^^E2 and using equation 
([8]) along with the equations satisfied by the background quantities e and Ti. 
gives an integral where the the integrand is a sum of terms each of which has 
a factor $2, |<i>'P, $$' or |V<i>p. The aim is to show that the sum of these 
terms is non-negative. To do this it is first assumed that the coefficient of |V$p 
is non-negative. This leads to the condition (P^j) . Next it is shown that the 
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quadratic form in $ and $' is positive semidefinite. This can be done by using 
the inequahty 

|AH$$'| < -A2h2$2 ^ 1 ($')2^ (34) 

2 2(5 

which holds for any (5 > 0, to estimate the quadratic form from befow by the 
foUowing sum of a term containing $^ and one containing 



2A 



1 



3/(e) 



2A- 6(1 + /'(£)) 



- 3 ( /'(e) 

A-i 



/(e) 



$2 



A" 



3 f{e) 



fie] 



1$' 



(35) 



It remains to ensure that the coefficients of these terms are non-negative and 
this foUows from ([30|) and (f3T|) , choosing S sufficiently small. It can be concluded 
from (j33|l that £2(1]) = 0(0(77)^^'^) as 77 ^ 0. As in the case of a linear equa- 
tion of state, corresponding estimates hold for spatial derivatives and pointwise 
estimates follow by Sobolev embedding. An integral formula for $' can be ob- 
tained as in the case of a linear equation of state. It reads (with some arguments 
suppressed; recall e and Ti. depend on rj) 



<^>'{V) = (/(e) + e) 



/(eo) + eo 
■1 I 



/(e) 



/'(e)A$ + 3 /'(e 



/(e) 



(36) 



If no further assumptions are made on the equation of state then using the 
known boundedness statements and repeatedly substituting into the right hand 
side of (j36p would lead to unwieldy expressions involving iterated integrals. 
Simpler results can be obtained if it is assumed that in the limit e — + 00 the 
function / is linear in leading order with lower powers as corrections. In other 
words for this assume that / admits an asymptotic expansion of the form 



fie) 



00 

E 



(37) 



as e — > 00. Here the fj are constants while {uj} is a decreasing sequence of real 
numbers all of which are less than one and which tend to —00 as j —> 00. Assume 
further that the relation obtained by differentiating this expansion term by term 
any number of times is also a valid asymptotic expansion. To have a concrete 
example, consider the polytropic equation of state which is given parametrically 
by the relations 

e = TO + Knm~^ , p — Km~^ (38) 

with constants K and n satisfying < K < 1 and n > 1. In this case the 
asymptotic expansion is of the form 



/(e) = n-ie 



+ n ^(ifn)"+ie"+i -|- 



(39) 
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Returning to the more general case ([37]), define a quantity m by 



m(e)=exp|^ + (40) 

Substituting the asymptotic expression ([57]) into (PU]) gives a corresponding 
asymptotic expansion for the function m(e) as a sum of powers of e with the 
leading term being proportional to e™+i . It follows from the continuity equation 
([U for the fluid that m is proportional to a^^. This leads to an asymptotic 
expansion for e in terms of a. The equation (O implies that a' = ^ySnG/Se^/^a^; 
substituting for e in terms of a gives rise to a relation which can be integrated to 
give an asymptotic expansion for a in terms of rj in the limit ?7 — > 0. The leading 

2 

term is proportional to ry3"'+i . Substituting this back in leads to an asymptotic 
expansion for a' from which an asymptotic expansion for Ti. can be obtained. An 
asymptotic expression for e in terms of 77 can also be derived. Thus in the end 
there are expansions for all the important quantities in the background solution 
in terms of rj. In all cases the leading term in the expansion agrees with that in 
the case of a linear equation of state. The result is an integral equation which 
can be written in the form 



C- h^iO^^iO -h:i{Q)HQ)dQ 
Jo 







(41) 



where C is a constant depending only on the data at time 770 and asymptotic 
expansions are available for the functions hi, /12 and h^. The leading terms in 
hi and /12 are constant multiples of the corresponding powers of r] for a linear 
equation of state. To see the leading order behaviour of /13 recall that in 
the coefficient of $ is 



(/'(e) - 



/(e) 



H. (42) 



Hence if fj is the first non- vanishing coefficient in the expansion (|39|) then the 

leading order power in h^ is less than that in by a = 2 — ^^"'""'"j^^g^""^'' . To 
obtain estimates close to rj = the estimate for the energy can be applied 
starting from 77 very small. In other words, cq can be chosen as large as desired. 
Then all the coefficients in the left hand side of ([50]) not involving A are as close 
as desired to those for the corresponding linear equation of state. Since A is 
arbitrarily small, the coefficient involving A is also arbitrarily small. It follows 
that A can be chosen to have any value strictly greater than 3(1 + w). Hence 
E can be bounded by any power greater than the power in the corresponding 
linear case. This is enough to proceed as in the proof of Theorem 1 to obtain 
an asymptotic expansion for $ where each invidual term is a constant multiple 
of an expression of the form 77'"'(log77)' with I = ot I — 1 and the leading 
term is just as in Theorem 1 with the corresponding value of w. The key thing 
that makes this work is that a < 2 so that no logarithms are generated when 
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evaluating the integral in ([H]) in the course of the iteration. The results of this 
discussion can be summed up as follows. 

Theorem 3. Let ^ be a smooth solution of ([5]) on (0, oo) x . Suppose that 
the equation of state has an asymptotic expansion of the form p7p . Then there 
are coefficients ^k,i with k > —2v belonging to an increasing sequence of real 
numbers tending to infinity and I € {0, 1}, smooth functions on , such that 
the formal series $fc^; (log 77)^77'^ is asymptotic to $ in the limit ry — s- 0. All 
coefficients in the expansion are determined uniquely by $-2i/,o o,iT-d ^0,0- 

4 Late-time asymptotics for a linear equation of 
state 

In this section information is obtained about the asymptotics of solutions of 
equation ([T]) in the limit rj —^ 00; some extensions of these results to more 
general equations of state are derived in Sect. [5l Once again energy estimates 
play a fundamental role. In this case it is convenient to treat homogeneous 
solutions separately. By a homogeneous solution we mean one which does not 
depend on the spatial coordinates. These can be characterized as the solutions 
whose initial data on a given spacelike hypersurface do not depend on the spatial 
coordinates. For this class of solutions equation ^ can be solved explicitly with 
the result that $ = A + Bi]^^^ for constants A and B. A general solution can 
be written as the sum of a homogeneous solution and a solution such that $ 
has zero mean on any hypersurface of constant conformal time. Call solutions 
of the latter type zero-mean solutions. Then in order to determine the late-time 
asymptotics for general solutions it suffices to do so for zero-mean solutions. In 
this case define ip{ri) — if^^^irj). Then ip satisfies the equation 

V'" = wAV^ +{^^ -''ij ^^^^ (43) 

Define an energy by 

E3{v) = l I Wm^ + M'^i^iv)?- (44) 

Then 

E'M = 2 _ 1^ ^-2 ^(^)^'(^). (45) 

The integral on the right hand side of this equation can be bounded, using the 
Cauchy-Schwarz inequality, in terms of the norms of ip' and ip. The first of 
these can be bounded in terms of the energy and due to the fact that the mean 
value of ip is zero, the same is true of the second. Thus E'^{ri) < Crj^^E^{rj) for 
a constant C. By Gronwall's inequality it follows that is globally bounded 
in the future. These arguments apply equally well to spatial derivatives of ip 
of any order. By the Sobolev embedding theorem it can be concluded that -0 
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and its spatial derivatives of any order are bounded. The energy bounds and 
the basic equation then imply that all spacetime derivatives of any order are 
uniformly bounded in time. 

Let rjj be a sequence of times tending to infinity and consider the trans- 
lates defined by ipj{r]) = tp{r] + ijj). The sequence Vj satisfies uniform C°° 
bounds. Consider the restriction of this sequence to an interval [f?o,r?i]. By the 
Arzela-Ascoli theorem the sequence of restrictions has a uniformly convergent 
subsequence. By passing to further subsequences and diagonalization it can be 
shown that ip and its spacetime derivatives of all orders converge uniformly on 
compact subsets to a limit W. Passing to the limit in the evolution equation 
for tjj along one of these sequences shows that W satisfies the flat-space wave 
equation W" — wAW. Note that a priori the function W could depend on the 
sequence of times chosen. This issue is examined more closely below. 

Given a smooth solution of ((43|) it is possible to do a Fourier transform in 
space to get the equation 

i/," = -w\k\^4> + (^1/2 - 1^ ,y-2^ (46) 

which is referred to below as the mode equation. Here /c is a vector. The 
restriction to zero-mean solutions implies that the case fc = of can be 
ignored. 

Lemma 1. Any solution (f> of equation (I46p has an asymptotic expansion of the 
form 

$irj) = WkC0siV^\k\{r^-fjk)) + 0i7j-^), (47) 
for constants fjk and Wk, in the limit rj oo. 

Proof. To prove the lemma it is convenient to introduce polar coordinates asso- 
ciated to the variables ijj and V''- Thus — rcos9 and --^=^^tp' = rsinO. 
This leads to the equations: 

1 



— , , , 1 - - rr)-^ sin cos (48) 
w|fc| V 4/ 

0' ^ -V^N + ^(^^-i).-^cos^^ (49) 

It follows from dH]) that 

0(77) = -V^|fc|(77-77fc) + 0(7,-1) (50) 
for a constant fjk- From (|48| it follows that 

r-(77) =VKfc(l + 0(77-1)) (51) 
for a constant Wk- As a consequence of ([50]) we have 

cos(7,(6l)) = cos(V^|/c| (7y - fjk)) + 0(77-^). (52) 
Together with ([5T|) this gives the conclusion of the lemma. □ 
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Consider a zero-mean solution of the type considered before. Let a function 
W be defined by taking the sequence rjj used above to consist of integer multiples 
of 27r. We now show that the function ip — W tends to zero as — > oo. In order 
to do this it suffices to show that it does so along a subsequence of an arbitrary 
sequence of values Q of rj tending to infinity. By passing to a subsequence 
as before it can be arranged that the translates by the amounts Cj converge 
uniformly on compact subsets as j — > oo. Call the limit Y. The aim is to 
prove that Y — 0. If not there must be some mode Y which is non-zero. It 
can be obtained as the limit of some '0 — W. From Lemma [T] it can be seen 
that W = Wk cos{^/w\k\{r] — rjk))- Hence ■>!) ~ W — 0{r]^^) and so F = 0, a 
contradiction. Convergence of derivatives can be obtained in a corresponding 
way. Thus any solution can be written as ^{r],x) = rj^''^'^ (^W{r], x) + o(l)). A 
similar result for the polarized Gowdy equation with a sharper estimate on the 
error term was proved in j4j. 

A late-time asymptotic expansion has now been derived which involves a 
solution W of the flat-space wave equation. Comparing with the results on 
parametrizing solutions by the coefficients in an asymptotic expansion near the 
singularity it is natural to ask if the function W can be prescribed freely. It will 
now be shown that this is the case by following the proof of an analogous result 
for the polarized Gowdy equation due to Ringstrom [11]. Write an arbitrary 
zero-mean solution in the form 

^{ri,x) ^ T]^''^^W{r],x) + uj{j],x). (53) 

Then uj satisfies the equation 

w" + i]~'^uj' - wAuj = ^1/2 - 1^ rj-^-^W. (54) 

Define 

H{V)^1 [ W{v)\'+w\WLu{fj)\' (55) 

and 

r(^) f coivWir,). (56) 

The aim is to study late times and attention will be restricted to the region where 
r] > w~^. At this point it is necessary to assume that w > 0. The following 
inequalities show the equivalence of H and iJ -I- F as norms of {uj' , Vw): 

TOI < ^^('7) , Ih<H + T<^H. (57) 

Now 



-f[H + rl = -i(H + r,-i^-±lr 
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Using the equivalence of H + T and H this can be used to derive the following 
differential inequality 




(H + T) 



1/2 



(59) 



(60) 



This quantity satisfies an inequality of the form 



(61) 



for a positive constant C depending on w. Since ri~^~^ is integrable at infinity 
this inequality can be used in just the same way as the corresponding inequality 
in [TT]. In this way it can be proved that given a solution W of the flat space 
wave equation there is a corresponding solution $ of ([1]). It follows from the 
proof that £4(7]) = 0(77"^""^). Hence H{r]) = 0(1]^'^"^^) and the solution 
decays like ry"'^/^^". 

The information obtained concerning the asymptotics of the solutions con- 
structed starting from a solution W of the wave equation is stronger that what 
was proved about general solutions of ([IJ up to this point. This can be improved 
on as follows. Given a solution $ of ([ij a solution W of the flat space wave 
equation is obtained. From there a solution $ of ([T]) is obtained with stronger 
information on the asymptotics. The aim is now to show that $ = <!>. To do 
this it is enough to show that each Fourier mode agrees. This means showing 
that a solution tp of (jlS)) vanishes if it tends to zero as 77 ^ 00. That the latter 
statement holds follows easily from (|48|) . What has been proved can be summed 
up in the following theorem. 

Theorem 4. Let ^ be a global smooth solution of ([1]). Then there exist con- 
stants A and B and a smooth solution W of the equation W" = wAW with zero 
spatial average such that 



This asymptotic expansion may be differentiated term by term in space as often 
as desired. 

Note that the third explicit term in this asymptotic expansion is often no larger 
than the error term. The function W can be prescribed freely. 

5 Late-time asymptotics for a general equation 
of state 

It will now be investigated how the results of the previous section can be ex- 
tended to the case of a more general equation of state. The class of equations of 



(62) 
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state which will be treated is defined by requiring that they admit an asymptotic 
expansion of the form 



oo 

f{e)^we + J2fj^''' (63) 



for e ^ 0. Here ui > 0, the coefficients aj are all greater than one and form 
an increasing sequence. To ensure the positivity of /' it is assumed that if 
w = the coefficient /i is positive. This form of the equation of state may be 
compared with that of ([57]) . It is further assumed that this expansion retains 
its validity when differentiated term by term as often as desired. An example 
is given by the polytropic equation of state (|55|) . In that case w — 0, fi — K 
and oi = ■^TT"- With this assumption information can be obtained on the 
leading order asymptotics of the background solution as 77 — *■ cxd. To simplify 
the notation define a — ai — I. It is convenient to use the mass density once 
more, writing (|40|) in the equivalent form 

m(e) = exp|-^\e + /(0)"'£^e| (64) 

Then m(e) has an expansion about e = where the leading term is proportional 
to e™+i . In particular, when w ^ the leading term is linear. Using the fact 
that TO is proportional to a^"^ for any equation of state leads to an asymptotic 
expansion for e in terms of a. Putting this information into ^ shows that 0(77) 

has an expansion in the limit 77 — )■ 00 with the leading term proportional to 

2 

7y3m+i . Finally it follows that e and Ti have expansions with leading terms pro- 

6(1 + ™) _ 

portional to rj 1+3™ and 77 respectively. With the leading asymptotics of the 
background solution having been determined it is possible to derive asymptotics 
for the coefficients in the equation for $. 

As in the case of a linear equation of state it is convenient to treat homoge- 
neous and zero-mean solutions separately. The homogeneous solutions will be 
analysed first. This leads to consideration of the equation obtained from ([8]) 
by omitting the term containing spatial derivatives. It is convenient here to 
exclude the case of a linear equation of state which was previously analysed so 
as to ensure that a is defined uniquely in terms of the equation of state. The 
coefficients satisfy: 

3(1 + /'(e)) = 3{l + w + {a+ l).hv~^) + o{r^) (65) 

and 

3 (/'(e) - ^) = 3/1'^'?-'' + o{r^) (66) 



where (3 — ^'^.^^^ . Define 



F = + a77-2-^$2 (g7) 
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where a is a positive constant which needs to be chosen appropriately in what 
follows. Computing the derivative of F with respect to rj and using the equation 
gives a sum of terms involving and The aim is to show that F is 

bounded and to do this it suffices to consider arbitrarily late times. The leading 
order terms in the coefficients of and are — ^rrr^r?"^ and —A'n~^~/^ 
respectively, where A is a positive constant. The coefficient of $$' has a leading 
term proportional to 77"^"'' for a general choice of a. However if a is chosen to 
be half the coefficient of the leading order term in the expansion of the coefficient 
of <f> in ([8]) then a cancellation occurs and the coefficient becomes o{t]~^~^). This 
choice is made here. The aim is to show that the term containing can be 
absorbed by the sum of the other two so as to leave a non-positive remainder. 
To do this the inequality 

is used. The powers of rj which arise from this inequality match those in the 
leading order terms in the coefficients of the manifestly negative terms in the 
expression for the derivative of F with respect to ij. Thus at late times the cross- 
term can be absorbed in the terms with the desired sign. The conclusion is that 
F is bounded. In fact this can be improved somewhat. The derivative of F can 
be estimated above by —2'^rj^^F for any positive constant 7 < + 1. This 
means that $' decays like ri~^ . It can be concluded that is bounded. From 
the evolution equation for $ and the boundedness statements already obtained 
it follows that (77^"^+^$')' is integrable. Thus $ = A + Brj'^" + . . . for constants 
A and B and the leading order behaviour is as in the case of a linear equation 
of state. 

It turns out to be useful for the analysis of the zero-mean solutions in the 
expanding direction to introduce a new time variable t satisfying the relation 
dr/di] = /'(e). Substituting the asymptotics of /'(e) in terms of rj into this 
provides an asymptotic expansion for r in terms of rj. For w > a linear relation 
is obtained in leading order while for w — and cr ^ ^ the expansion reads 

T = Cl 77^-3"+ Too + ... (69) 

for constants Ci and Tcc . Note that the second term in this expansion is only 
smaller than the first for cr < i. For w — and a — ^ the power in this 
expression gets replaced by logr/. From these facts it can be seen that r — > cx) 
for rj ^ 00 when w > or when w = and cr < ^. In contrast r tends to the 
finite limit Too for — s- 00 when w = and a > ^. This is a symptom of a 
bifurcation where the asymptotics of the linearized solution undergoes a major 
change. For convenience we say that the dynamics for an equation of state with 
an asymptotic expansion of the form ([55]) is underdamped if w > or cr < i, 
critical if w = and f = i and overdamped if w = and a > ^. 

Next the late-time behaviour will be analysed for zero-mean solutions with 
an equation of state corresponding to underdamped dynamics. The first step is 
to introduce the time variable r into ((8]) with the result: 

+ 3Zn<i>r + 3Y'H^<P - A* = (70) 
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where 



r = f'{,)^M, (71) 

and Ti. — a^^Ur- Derivatives with respect to r are denoted by subscripts. Next 
the term containing <I>t- will be eliminated by multiplying <i> by a suitable factor 
fl~^. Choose rt to satisfy 

For all three types the behaviour of as a function of a in the limit e ^ can 
be determined. The result is that the leading order term in il is proportional 
to 0^5(1+"') for w > and proportional to a^2(i^f ) for w = 0. The function 
^ = ri~^$ satisfies an equation of the form 

= A{e)'HH + A* (74) 

where A{e) is a rational function of e, /(e), /"(e) and /'"(e). Under the given 
assumptions on the equation of state it is bounded. Proving this requires ex- 
amining many terms but is routine. For example the only term containing the 
third derivative of / is '^^^"'"^/'(e/ — ~- leading order terms in the asymp- 
totic expansions of numerator and denominator are both proportional to e*^. 
Note also that the leading order term in the expansion for Ti. is proportional to 
for any cr < i. Note for comparison that H tends to a constant value as 
T — * cxD in the case a — ^. 
Define an energy by 

E,ir) = ^l^l + \V^\^. (75) 

Then using the same techniques as in previous energy estimates shows that 
there is a constant C such that 

^ < ClAln^E, (76) 

UT 

Using the information available concerning A and Ti. shows that is bounded 
in the future. Taking derivatives of the equation and using the same arguments 
as in previous cases shows that 5* and its derivatives of all orders with respect 
to X and r are bounded. It follows that any sequence of translates 5'(t + t„) for 
a sequence t„ tending to infinity has a subsequence which converges on compact 
subsets to a limit W. 

Doing a Fourier transform of the equation ([74]) in space leads to the mode 
equation 

^r^^ = -|fc|2^-h ATY^* (77) 
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Introducing polar coordinates in the ( ^)-plane leads to the system 
dr 1 - 

_ = —ArU-'^ sine cos 9 (78) 
ctr |fc| 

^ = + J-AW-^cos^^ (79) 

dr \k\ 



This implies that 
and 



e{T)=-\k\{T-fk)+0{T-^) (80) 

r{T)=Wk{l + 0{T-^)) (81) 



for some constants fk and Wk- Arguing as in the case of a Hnear equation of 
state leads to the relation \l/(r, a;) = W{t,x) + o(l) where is a solution of 
the equation W^r = AW^. Using the form of the leading order term in f2 as a 
function of a, it can be shown that the leading order term in O as a function 

_Mi±wl 1 3(1- f) 

of T is given by r 1+^™ = t 2 and t "n~^W in the cases w > and w = 
respectively. Note that the first of these reproduces the result in the case of a 
linear equation of state. It does not seem to be possible to write the expansion 
directly in terms of r] in such a way that it gives more insight than the expression 
in terms of r. The leading asymptotics of a zero-mean solution is obtained by 
taking a solution of the flat space wave equation, distorting the time variable 
by a diffcomorphism and multiplying by a power of the time coordinate which 
has been explicitly computed. 

Consider next the case w = 0, a > ^ (overdamped case). The time coor- 
dinate T tends to the finite limit Tqc as ry ^ 00. Define G = THT^dr'H- The 
function G tends to the limit |(c7 — ^) as r ^ Too- Let 

Eq= y $2 + |V$|2 + A2w2$2 (82) 

for a constant A which remains to be chosen. For a constant A computing 
dria^^Ee) gives rise to a sum of expressions containing $2^ |V$p. 
Using the inequality 

n'^^r^ < Trfi^l + ttH^^^ (83) 
2o 2 

leads to an inequality where the term involving has been eliminated. To 
obtain some control on the energy by means of the inequality the coefficients A 
and A should be chosen in such a way that all terms on the right hand side are 
manifestly non-positive. The conditions for this to happen are the inequalities 
A < 0, 

-'-lA^-SFI <3Z-A (84) 
2o 

and 

^|A2-3y| < -A2(A + G) (85) 
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Note that these inequahties imply in particular that A < 0. Consider now the 
limit T ^ Too where Y behaves asymptotically like /icre°" and Z ^ 1 — ^a. In 
this limit the inequality ([85)1 reduces toA<— |(cr— -j) — |. Suppose therefore 
that A < — |(ct— i). Then by choosing S sufBciently small it can be arranged 
that the limiting inequality is satisfied. In the limit the inequality (|84p reduces to 
^ < 3— 1(7 — A. Choose A so that this inequality is satisfied strictly. With these 
choices both inequalities are satisfied strictly in the limit. For r sufhciently close 
to Too the coefficients in (|84)) and (|85p are as close as desired to their limiting 
values. Making them close enough ensures that these two inequalities continue 
to be satisfied. It follows that with these choices of the parameters dria^^EQ) 
is non-positive at late times. It can be concluded that = 0{a^'^^). This 
gives a limit on the growth rate of Eq in terms of that of the scale factor. As in 
previous cases corresponding estimates can be obtained for derivatives and as a 
consequence pointwise estimates derived. It follows that $ = 0{a~^Ti.~^). From 
what is known about the background solution it follows that TL is proportional 
to a5('^^ 3). Thus if p = — |(o- — i) — A then $ = 0{aP). This power is positive 
but may be made as small as desired by choosing A suitably. By the usual 
methods similar bounds can be obtained for spatial derivatives of 

To get more information about the asymptotics as r ^ Too it is convenient 
to rewrite the equation in terms of the new time variable s = Too — t. The 
resulting equation is 

^ss ' izn^s + 3r7i:2$ - A* = o (86) 

As s ^ the coefficient Z tends to 1 — |- while Ti. and YTi.^ are proportional 

_1 2 

in leading order to s and s^'''"^/^) respectively. The last exponent is positive 
for any a > 1/3 so that the corresponding coefficient tends to zero as s — > 0. 
Let B be a positive solution of the equation ^ = —iZJ-LB. Then ([5^ implies 
the following integral equation: 

$^ === 1 (^$1 +^ S(-3yH2$^^$)^ (87) 

for a function $i(x). Here the fact has been used that the integral occurring 
in this equation converges. This follows from the fact that in leading order 
B is proportional to s "-^/^ and the bounds already obtained for <I> and its 
derivatives. When B^^ diverges faster than s^^ in the limit s 0, which 
happens for cr < |, the known bounds on <& imply that $i =0. Hence $ is 
bounded in the limit s ^ in that case. When cr > | it can also be concluded 
that $ is bounded. For a — a logarithmic divergence of is not ruled out. In 
all cases the integral equation can be used to obtain an asymptotic expansion 
for <&. Schematically this expansion is of the form 

$(,7,a:) =^<l),(x)C,(?7). (88) 

i 

for some functions with Ci+i{v) — for each i. This is very different 

from the expansion in the limit ij ^ oo obtained when it; > or cr < i. In the 
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present case, scaling the solution by a suitable function of r] gives a result which 
converges to a function of x as 77 — > cx). In the other case a similar rescaling 
can lead to a profile which moves around the torus with constant velocity. (In 
general it leads to a superposition of profiles of this kind.) In the latter case 
there are waves which continue to propagate at arbitrarily late times. In the 
case c > 5 the waves 'freeze'. This is reminiscent of the late-time asymptotics 
of the gravitational field in spacetimes with positive cosmological constant (cf. 
[TO]). 

To make this argument more concrete consider the special case where the 
equation of state is /(e) = /ie°'"'"^ for some tr between ^ and |. Using the 
convergence of the integral it follows immediately that $(s, x) = ^o{x) + O(s^) 
for some function $9. Putting this information back into the integral equation 
gives $(s,x) = $o(a;) + ^^A<i>o{x)s^ + .... 

Consider finally the case w = 0, a = ^ (critical case). Then rj = rjoe'^ + . . . 
where 770 is a constant and Ci corresponds to the constant appearing in (|69p . 
The arguments leading to the estimate $ — 0{tP) can be carried out as in 
the case ct > i. The only difference is that the limit t ^ t^o is replaced by 
T ^ 00. In the case cr = i the quantity tends to a constant for t — > 00 

and YH'^ is proportional to e^'^^ in leading order. A quantity B can be 
introduced as before and an integral equation obtained. In this case _B is a 
decaying exponential. Unfortunately it does not seem to be possible to use this 
integral equation to refine the asymptotics in this case and this matter will not 
be pursued further here. 

Acknowledgements 

The authors gratefully acknowledge the hospitality and financial support of the 
Mittag-Leffler Institute where part of this work was carried out. 

References 

[1] Andersson, L. and Rendall, A. D. 2001 Quiescent cosmological singularities. 
Commun. Math. Phys. 218, 479-511. 

[2] Chandrasekhar, S. 1961 Hydrodynamic and hydromagnetic stability. Ox- 
ford University Press, Oxford. 

[3] Isenberg, J. and Moncrief, V. 1990 Asymptotic behaviour of the gravita- 
tional field and the nature of singularities in Gowdy spacetimes. Ann. Phys. 
199, 84-122. 

[4] Jurke, T. 2003 On the future asymptotics of polarized Gowdy T'^-models. 
Class. Quantum Grav. 20, 173-191. 

[5] Keller, E. F. and Segel, L. A. 1970 Initiation of slime mold aggregation 
viewed as an instability. J. Theor. Biol. 26, 399-415. 



23 



[6] Kichcnassamy, S. 1996 Fuchsian equations in Sobolev spaces and blow-up. 
J. Diff. Eq. 125, 299-327. 

[7] Kichenassamy, S. and Rendall, A. D. 1998 Analytic description of singu- 
larities in Gowdy spacetimes. Class. Quantum Grav. 15, 1339-1355. 

[8] Mukhanov, V. 2005 Physical foundations of cosmology. Cambridge Univer- 
sity Press, Cambridge. 

[9] Rendall, A. D. 2000 Fuchsian analysis of singularities in Gowdy spacetimes 
beyond analyticity. Class. Quantum Grav. 17, 3305-3316. 

[10] Kendall, A. D. 2004 Asymptotics of solutions of the Einstein equations with 
positive cosmological constant. Ann. Inst. H. Poincare 5, 1041 1064. 

[11] Ringstrom, H. 2005 Data at the moment of infinite expansion for polarized 
Gowdy. Class. Quantum Grav. 22, 1647 1653. 

[12] Turing, A. M. 1952 The chemical basis of morphogenesis. Phil. Trans. R. 
Soc. Lond. B237, 37-72. 



24 



